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there  is  the  "diabatic"  potential  surface  formed  from  the 
electronically-uncoupled  reactants  and  environment  and  the 
diabatic  surface  formed  from  the  electronically-uncoupled 
products  and  environment.  The  electronic  interaction  between 
A and  B introduces  new  surfaces — the  adiabatic  potential  energy 
surfaces,  discussed  for  example  by  the  principal  investigator 
in  Ann.  Rev.  Phys.  Chem. , 15^  155  (1964).  It  is  desirable  to 
develop  a numerical  method  to  treat  the  dynamics  of  molecular 
motion  on  those  surfaces,  to  calculate  the  reaction  rate.  In 
the  present  Technical  Report  Dr.  Babamov  formulates  such  a 
method.  In  practice,  to  reduce  the  number  of  quantum  states 
needed  in  the  computation  one  would  need  to  simplify  the  actual 
many-dimensional  problem  before  applying  the  method.  One  such 
application  has  been  made  by  Drs.  Babamov  and  DiGiacomo  of 
this  laboratory. 
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ABSTRACT 


A new  method  for  solving  the  coupled  equations  in  the  adiabatic 
representation  which  arise  in  quantum  mechanical  treatment  of  atomic  and 
molecular  collisions  Is  developed.  The  solutions  of  the  uncoupled 
adiabatic  equations,  obtained  from  the  coupled  equations  in  the  adiabatic 
representation  by  neglecting  the  coupling  terms,  are  used  to  derive  a set 
of  new  exact  equations  for  the  purely  Inelastic  effects  of  the  collision. 
The  solutions  of  the  new  equations  are  practically  nonoscl I latory  and 
can  be  solved  numerically  with  little  computational  effort. 
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I.  INTRODUCTION 

The  computational  wave  mechanical  treatment  of  inelastic  atomic  and  mole- 
cular collisions  is  frequently  reduced  to  solving  a set  of  coupled  ordinary 
linear  differential  equations.'"5  The  coupled  ordinary  differential  equations 
••  are  obtained  by  expanding  the  solution  of  the  t ime- independent  Schrodinger 
equation  for  the  system  in  a product  basis  set  consisting  of  known  wavef unc- 
tions for  the  internal  motion  and  for  the  angular  relative  motion,  multiplied 
by  unknown  coefficients  which  depend  on  the  relative  radial  motion  coordinate. 

The  choice  of  the  basis  set  functions  for  the  internal  motion  determines 
the  structure  of  the  resulting  coupled  equations.  When  the  wavef unct ions  for 
the  internal  motion  of  the  separated  particles  are  used  in  the  basis  set  the 
coupled  equations  are  said  to  be  in  the  diabatic  representation.  On  the 
other  hand  when  the  solutions  of  the  Schrodinger  equation  for  the  colliding 
pair  at  each  fixed  interparticle  distance  are  used  in  the  expansion,  the  resulting 
equations  are  In  the  adiabatic  representation.3-6 

The  coupled  equations  in  the  diabatic  representation  are  used  often 
In  the  treatment  of  vibrat ional 1 y and  rotationally  inelastic  molecular 
collisions  occurring  on  a single  electronic  potential  energy  surface.7** 

Developing  methods  for  solving  the  coupled  equations  in  the  diabatic  repre- 
sentation has  thereby  received  considerable  attention,  particularly  In  the 
past  decade,  and  a number  of  highly  efficient  methods  for  solving  them  have 
been  developed. 7-1 0 

The  equations  In  the  adiabatic  representation  are  frequently 
encountered  when  a rigorous  treatment  of  an  electronically  nonadiabatic  atomic 
collision  is  attempted.  ' ' *' * They  have  a somewhat  more  complex  mathematical 
structure  and  developing  methods  for  solving  them  has  received  relatively 
little  attention  In  the  I Iterature.1  3 * 14  The  most  widely  used  of  the  new7"* 


efficient  methods  for  solving  the  equations  in  the  diabatic  represen- 
tation hav*e  not  been  adapted  for  solving  the  adiabatic  equations.  Although 

a direct  numerical  integration  of  the  equations  in  the  adiabatic  representation 
Is  possible  and  has  been  conducted  for  some  model  probl ems 1 3* 1 4 the  only  method 
used  has  been  direct  numerical  integration  of  the  coupled  equations  based 
on  a polynomial  approximation  to  the  wavef unct ion.  The  inelastic  atomic 
collisions  have  in  practice  usually  been  treated  instead  by  solving  a set  of 

diabatic  equations  constructed  in  some  approximate  manner,15  or  still  more 

approximately  by  solving  the  corresponding  so-called  classical  path  equations.  1 

The  latter  equations  are  obtained  by  treating  the  electronic  motion  quantum 

mechanically  and  the  relative  motion  of  the  atoms  classically. 

In  this  paper  a stable  and  efficient  method  for  solving  a set  of  coupled 
ordinary  differential  equations  in  the  adiabatic  representation  is  outlined. 

The  method  can  be  used  for  treating  electronic  transitions  and  charge  trans- 
fer in  atomic  collisions,  which  can  be  rigorously  formulated  in  the  adia- 
batic representat ion5* 1 * The  electronic  transitions  in  molecular  collisions 
can  also  be  treated  by  the  present  method  since  they  are  also  conveniently 
formulated  in  the  adiabatic  representation.  13 

Both  the  spatial  extension  of  the  coupling  region17  and  the  number  of 
coupled  states'*  are  smaller  in  the  adiabatic  .representation  than  in  the 
diabatic  one.  Since  the  equations  in  the  diabatic  representation  and  the 
equations  in  any  other  representation  in  which  both  diabatic  and  adiabatic 
coupling  elements  appear  can  be  straightforwardly  converted  into  a pure 
adiabatic  representation3.*5  the  present  method  could  also  be  used  io  cases 
when  a diabatic  -or  some  other  representation  is  most  convenient  for  Initial 
formulation  of  the  problem 
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The  Hamiltonian  for  the  system  and  the  usual  quantum  mechanical  coupled 
equations  are  given  in  Sec.  II.  The  structure  and  the  symmetry  properties  of 
the  coupled  equations  are  analyzed  in  the  Appendix. 

“ The  solution  of  the  coupled  equations  is  performed  in  three  stages.  First, 
a set  of  uncoupled  second  order  differential  equations,  obtained  from  the 
equations  in  the  adiabatic  representation  by  neglecting  the  coupling  matrix 
elements,  is  solved  using  one  of  the  standard  methods  in  the  literature.* 

The  solutions  of  these  uncoupled  equations  are  then  used  in  Sec.  Ill  to 
convert  the  equations  In  the  adiabatic  representation  into  a set  of  new 
first-order  coupled  differential  equations  by  means  of  a variation  of  con- 
stant type  of  treatment.  The  dependent  variables  in  these  new  equations  are 
typically  much  more  slowly  varying  than  the  original  ones.  The  resulting 
first-order  equations  (Eq.  3.14)  can  then  be  solved  using  one  of  the  standard 
methods  to  evaluate  the  scattering  matrix  elements. 

The  method  is  illustrated  in  Sec.  IV  using  as  an  example  the  evalu- 
ation of  inelastic  transition  probabilities  for  a He+  + Ne  collision  in 
the  two-state  approximation.  The  model  problem  used  has  been  treated 
previsouly  by  other  methods,  and  so  the  feasibility  of  the  present  treat- 
ment can  be  verified. 


* 
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II.  THE  COUPLED  EQUATIONS 

The  Schrodinger  equation  for  a system  of  two  colliding  particles  which 
have  internal  structure  is  given  In  reduced  units  by: 

- E V*  + VU, R)]  Y=  ET  . (2-1) 

where  the  vector  R is  defined  by  its  magnitude  R and  the  two  angles  describing 

its  orientation  denoted  collectively  by  n ; R denotes  the  relative  position 

of  the  centers  of  mass  of  the  two  particles.  The  vector  £ denotes  all 

internal  coordinates.  The  sum  7* is  over  all  the  internal  coordinates. 

The  solution  of  the  time- independent  Schrodinger  equation  with  the 

Hamiltonian  (2.1)  can  be  expanded  in  a product  basis  set  consisting  of  a 

set  of  known  wavef unct ions  r • H ;R)  which  depend  on  r and  n for  every  value 

12d 

of  R,  multiplied  by  R-dependent  coefficients: 

4 

Y = E ♦,(R)  X,  (r,n;R)  (2,2) 

It  is  also  convenient  to  choose  X-  as  a set  of  real  functions. 

Substituting  (2.2)  Into  the  Schrodinger  equation  and  integrating  over  £ 
and  0 a *et  of  coupled  ordinary  differential  equations  for  the  radial 
coefficient  wavefunctlons  f,(R)  Is  obtained.  s.  *•  »* 

p'V*  [2WbuVvij*j1  * E*' 

where 

p « -1&/&R.  P,j  * <X||p|Xj>.  ®,j  • <Xj!p*IXj> 

VU  - VJI  - <x,l"Mxj> 


(2.3) 


(2.4a) 


(2.4b) 
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H*  contains  the  entire  potential  plus  the  kinetic  energy  terms  for 

the  internal  coordinates  r and  for  the  orientational  coordinates  Q.  If  a 

^ aw 

rotating  .coordinate  system  is  used  H'  also  contains  all  the  terms  steaming 

from  the  rotation  of  the  coordinate  system.  12  The  brackets  < ) denote 
integration  over  jq  and  q 

mm  MW 

Eq.  (2.3)  can  be  written  compactly  In  matrix  notation: 

[Ip2  + 2Pp  4 B + V]  if  = E$  , (2.  5) 


where  4 is  a column  vector  with  elements  tj 
matrices  with  elements  given  by  (2.4),  while 
Applying  p to  the  scalar  product  of  two 
equals  6.j  one  has 


and  P,  B,  and  V ere  square 
the  symbol  I denotes  unit  matrix. 
X,'s.  and  noting  that  <x5 |Xj> 


0=  P<X,!Xj>  = * (pX,Kj>  + (X.fpKj)  = -Pjj  + P,j  , (2.6) 


Thus,  P is  a hermitian  matrix.  If  the  basis  set  X;  is  complete  the  matrix  B 
can  be  expressed  in  terms  of  the  matrix  P and  its  derivative3 
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| = P2  + (pp),  (2.7) 

where  the  operator  p is  not  considered  to  act  beyond  the  close  of  the  paren- 
theses, and  the  elements  of  P2  are  obtained  from  those  of  P by  matrix  multi- 
plication. Using  (2i7)  the  coupled  equations  (2.5)  become 

(T  + V)  * = (2.8) 

T = Ip2  + 2Pp  + (pP)  + P*  = (Ip  + P)2  (2.9) 

It  is  In  practice  easier  to  evaluate  pP  and  P2  In  (2.9)  than  the  B in  (2.5)  and  (2.4) 
since  evaluation  of  | requires  numerical  evaluation  of  second  derivatives 
of  X,. 


Eq.  (2.  8)  is  a set  of  ordinary  differential  equations,  usually  an  Infinite 
set,  end  is  a partly  integrated  form  of  the  Schrodinger  equation,  fully 
equivalent  to  (2.1).  The  matrix  operator  on  the  l.h. s.  of  (2. B)  then  represents 
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the  Hamiltonian  operator  H from  (2.1)  as  applied  to  the  space  of  coefficient 
wavefunct ions  i(f.(R).  The  matrix  operator  T in  (2.9)  arises  from  the 
radial  part  of  the  kinetic  energy  operator  in  (2.1)  and  can  be  viewed  as 
the*  radial  kinetic  energy  part  of  the  matrix  Hamiltonian  operator. 

The  matrix  V represents  an  effective  matrix  potential  for  the  radial 
motion. 

It  is  convenient  for  simplifying  the  following  discussion  to  consider 
instead  of  Eq.  (2.8),  which  involves  the  vector  ^ the  following  matrix 
equation: 

[Ip*  + 2Pp  4 (pf)  4 P*  4 y]j  = E|  (2.  10 ) 

where  j|  denotes  an  n by  n square  matrix  which  contains  as  columns  n linearly 
independent  vector  solutions  4 of  (2.6). 

The  matrix  V In  (2.10)  can  be  rewritten  as 

V - V°  4 Vint  ; V°j  - 6i}  (E|  4 ijU^O/R2)  (2.11) 

where  e.  is  the  eigenvalue  of  the  internal  (^-dependent)  part  of  the  basis 
function  Xj » and  is  the  angular  momentum  associated  with  the  angular  (jg 
dependent)  part  of  the  Xj • In  the  limit  of  large  R the  adiabatic  coupling 
matrix  P and  the  potential  coupling  matrix  y*nt  vanish  and  the  coupled 
equations  (2.11)  can  be  written  as  a set  of  uncoupled  equations 

(jp2  + v°)ik  “ ; (2.12) 

the  diagonal  angular  terms  £j(£|4l)/R2,  which  also  vanish  for  large  R,have  been 
retained  in  (2.12)  to  simplify  the  structure  of  the  phases  in  the  scattering 
matrix.  The  solution  of  (2.12),  which  for  large  R is  equivalent  to  (2.11)  is 
given  in  terms  of  Riccati-Bessel  functions  and  can  be  written  for  large  R using 
the  assymptotic  form  of  the  Riccati-Bessel  functions  as: 


J(R)  - £*(R)£+  * F-(R)C",. 


(R  - -) 


. (2<13) 
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where  Ft  ere  diagonal  matrices  with  elements 

Fyj (R)  = ( l/2k<)  */*  exp[+i(k.R-i.Tr/2]  (2.  l4) 

l.e.  Fj*  and  Fjj  are  the  diagonal  elements  of  the  two  1 inearl y-independent 
;esymptotic  solutions  of  (2.  14).  The  factor  (1/2)'^*  has  been  inserted  to 

4.  • «•  % 4 . 

normalize  them  to  a Wronskian,  F. jFj .-Fj .F. j , equal  to  -i,  where  the  dot 
denotes  d/dR.  The  ci  in  (2.13;  are  integration  constant  matrices.  The 
choice  of  the  phases  in  (2.14)  is  a matter  of  convention;  different 
phases  will  be'used  in  Sec.  3 to  define  the  purely  inelastic  scattering 
matrix. 

The  scattering  matrix  (S-matrix)  is  usually  introduced  as  the  matrix 
that  converts  C”  in  (2.13)  into  C+I,i. 

?+  = S C'  (2-15) 

The  choice  of  the  basis  set  X.(r ,0)  in  (2.2)  determines  the  structure 

of  the  matrices  in  (2.10;.  Two  particular  choices  for  which  the  structure 

of  (2.10)  is  qualitatively  simplified  are  of  practical  interest.  If  the  basis 
set  is  chosen  to  be  the  set  of  eigenfunctions  of  the  Internal  orientational  part 
of  the  Hamiltonian  for. large  R the  matrix  P vanishes  identically  (due  to  the  lack 

of  dependence  of  the  basis  set  on  R)  and  all  the  coupling  In  the  differential 

equation  is  contained  in  the  effective  potential  energy  matrix  V.  Such  basis 
set  is  a dtabatlc  basis  set  and  the  resulting  set  of  coupled  equations  is 
In  the  dlabatic  representation.4 


The  coupled  equations  in  the  dlabatic  representation  are  of  the  form 

(IP*+V)J  «=  EJ  (2.16) 
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If  the  basts  set  ts  chosen  Instead  to  be  the  set  of  eigenfunctions  of  the 
internal -orientational  part  of  the  Hamiltonian,  i.e.,of  H1  at  every  R,  the  of f-d iagonal  i 
elements  of  V(2.4b)  vanish  identically,  sc  that  all  the  coupling  is  dynamic,  and 
lt-fs  contained  in  the  kinetic  energy  part  of  the  Hamiltonian.  The  diagonal 
elements  of  V then  obey  the  well  known  non-crossing  rule,19  i.e. , they  cannot 
cross  at  any  R unless  the  matrix  element  of  P that  couples  them  vanishes. 

In  the  adiabatic  representat ion,  the  equations  are  of  form: 


Ijp*  + 2Pp  + (pP)  + P*  + V ]|  = EJ 


where  V is  a diagonal  matrix  with  elements. 


Vij  " 6ij€i> 


(2.17) 


(2.18) 


and  Cj  is  the  i-th  eigenvalue  of  H'  for  a given  R.  All  the  coupling  is  contained 
in  the  kinetic  energy  matrix,  which  in  turn  is  given  in  terms  of  the  adiabatic 
coupl ing  matrix  P. 

It  should  be  noted  that  in  some  cases  (e.g.,  when  the  angular  coupling  is 
not  negligible)  the  eigenfunctions  of  the  internal-orientational  part  of  the 
Hamiltonian  H1  for  every  R may  be  difficult  to  evaluate  directly.  In  such 
cases  using  an  approximation  to  the  adiabatic  basis  set  (e.g.,  product  of 
eigenfunctions  of  the  internal  part  of  H multiplied  by  spherical  harmonics) 
would  lead  to  nonzero  of-diagonal  elements  in  ^ and  the  coupled  equations  (2.10) 
will  not  be  in  the  adiabatic  representation.  The  coupled  equations  can  then 
be  converted  to  the  adiabatic  representation  by  diagonalizing  the  potential 
energy  matrix  V at  every  R,  and  transforming  the  kinetic  energy  matrix 
accordingly.  The  same  can  be  done  when  the  diabatic  representation  can  be 
formulated  and  it  is  easier  to  set  up.  The  transformation  matrices  needed 
to  convert  the  equations  between  the  adiabatic,  diabatic  and  any  other  general 
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representation  are  given  in  the  appendix.  A rigorous  discussion  of  the 
symmetry  properties  of  the  coupled  equations  is  also  given  in  the  appendix. 

Once  a set  of  equations  in  the  adiabatic  representation  has  been  con- 
structed, either  directly  or  by  transforming  the  equations  in  some  other 
representation  the  method  described  in  the  next  section  can  be  used  to  con- 
struct their  solution. 
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III.  EQUATIONS  FOR  THE  AMPLITUDE  COEFFICIENTS 

Using  Eq.  (2.9),  Eq.  (2.10)  can  be  rewritten  as 

Qp  + P)*i  = F*i  , TT,J  = (E-VIJ)6Ij.  (3.1) 

For  most  8tomic  collisions  the  elements  of  the  coupling  radial  momentum 
matrix  P in  (2.10)  are  small  throughout  the  collision  range.  In  such 
cases  a good  approximation  to  the  solution  of  (3.1)  can  be  obtained 

by  neglecting  the  coupling  matrix  p and  solving  the  resulting  set  of  uncoupled 
second  order  differential  equations. 

IP*fe=n^  (3.2) 

Equation  (3.2)  can  be  solved  numerically  more  easily  than  (3.1); 
some  good  analytic  approximations  also  exist.2  0 Approximating  the  solu- 
tion of  (3.1)  by  a solution  of  (3.2)  Is  for  a system  of  two  atoms  essentially, 
the  Born-Oppenhe imer  approximation.11  It  corresponds  physically  to  the  neglect 
of  the  effect  of  the  relative  atomic  motion  on  the  electronic  states  of  the 
system.  For  slew  atomic  collisions  the  Born-Oppenhe imer  approximation  is 
usually  already  a very  good  one  except  for  one  important  special  case:  If 

t 

two  or  more  matrix  elements  of  the  diagonal  matrix  V are  nearly  degenerate 
for  some  values  of  the  radial  coordinate  R the  corresponding  matrix  elements 
In  the  matrix  P become  often  large  In  the  same  region  and  cannot  be  neglected 
any  more.  The  occurrence  of  such  avoided  crossing  of  the  diagonal  matrix 
elements  of  the  matrix  V as  a function  of  the  radial  coordinate  is  the  main 
reason  for  electronic  transitions  In  slow  atomic  collisions.  The  breakdown  of 
the  Born-Oppenheimer  approximation  in  this  case  is  localized  to  a relatively 
narrow  spatial  region  around  the  avoided  crossing;  Eg.  (3.2)  is  however, 
still  an  excellent  approximation  to  (3.1)  outside  the  avoided  crossing  region. 
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The  solution  of  (3.1)  subject  to  a given  initial  condition  at  small  R 
can  be  factored  at  any  point  into  the  solution  of  (3.2)  subject  to  the  same 
initial  condition  and  a set  of  additional  correction  factors  which  reflect 
the  nonvanishing  of  the  matrix  P. 

The  main  features  and  the  typical  oscillatory  structure  of  the  solution 
of  (3.1)  are  already  contained  in  the  solution  of  (3.2).  The  additional 
correction  factors  vary  more  slowly  with  R than  the  solution  of  (3.1)  itself 
and  are  practically  constant  outside  the  region  where  the  matrix  P is  large. 
In  the  remainder  of  this  section  the  differential  equations  for  the  above 
correction  factors  are  derived. 

Any  solution  of  each  of  the  uncoupled  equations 

p**il  = "*1  \t> j j (3*3) 

can  be  written  as  a linear  combination  of  two  1 Inearl y-independent  solutions 
of  (3.3)  denoted  by  ut.  and  U.., 

#ii  = UTiCTi  + UiiCii  ' (3*4) 

The  two  1 inearly  independent  sol  ut  ions  of  (3.  3) , lit.  and 

Uj j#  will  be  taken  to  be  the  solutions  with  the  asymptotic  form 

Ufi  ~ [1/(2^)'/*] exp  [±!(k,R  - 1 j)]  . (3.5) 

and 

which  are  complex  conjugates  of  each  otherAhave  been  normalized  in  such  a 
way  that  their  Wronskian,  Uf jUTj  - ujjU^j,  is  equal  to  -i. 

The  general  solution  of  (3.2)  can  then  be  written  In  form  of  a linear 

combination  of  the  two  solutions  U*  and  U". 

s c 

^ U+C++  u”c*  • (3.6) 


n 


The  only  nonvenishlng  elements  of  the  matrices  U are  the  diagonal 
ones  defined  by  (3.5),  while  the  c£  are  diagonal  constant  matrices. 

One  seeks  now  a solution  to  the  coupled  adiabatic  equations  (3.1) 
1 o terms  of  the  solutions  of  (3.2)  U+  and  U", 


| = uV(R)  + U"a"(R) 


(3.7) 


with  the  coefficients  a|j  being  functions  of  R.  Since  the  matrices  U“  are 
diagonal, the  elements  of  ^ are  given  by 


♦ jj  = ^iiaij^^  + unrt|j(R)> 


(3.8) 


Since  there  are  twice  as  many  ct^j  In  (3.7)  as  there  are  wave  func- 
tions fjj  from  which  the  ayj  can  be  determined,  one  is  free  to  impose  n* 
constraints  In  form  of  n*  algebraic  or  differential  relations  between  the 
a^j,  where  n Is  the  size  of  the  matrices  In  (3.1).  While  the  n*  constraints 
are  completely  arbitrary  at  this  point,  the  choice  of  the  constraint  will 
affect  the  ease  with  which  the  resulting  equations  will  be  solved.  The  most 
slowly  changing  coefficients  ay.  result  from  a constraint  that  makes  the  ay. 

J J 

come  closest  to  being  an  average  of  the  absolute  value  of  the  rapidly  oscil- 
lating wavefunction  solutions  of  the  coupled  equations.  In  the  classical 
limit,  when  the  period  of  oscillations  of  the  wavefunctions  goes  to  zero  the 
a^j  become  a measure  of  the  probability  of  finding  the  particle  in  the 
channel  1.  for  a system  initially  in  channel  J,  since  If  one  integrates  the 
absolute  square  of  the  wavefunction  in  a small  region  the  value  will  be  propor- 
tional to  the  value  of  the  absolute  value  square  of  the  much  smoother  functions 


Comparing  (3.1)  and  (3.2)  It  can  be  seen  that  applying  the  operator 
F = (ip+P)  twice  to  the  solution  of  (3.1)  has  the  same  effect  as  applying  the 
operator  p twice  to  the  solutions  of  (3.2).  The  operator  F*# 
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F * - tlP-P)*  . (5.9) 

when  applied  on  the  solutions  of  (3.2)  thus  acts  effectively  only  on  the 
elastic  part  of  the  solution,  given  by  the  solutions  of  (3.3  ),  U*,  but  not 
on  the  additional  factors,  rj*. 

f2u+^+  + y*g_  = (p2y+)s+  + (p2u")g-  (3.10) 

Defining  the  operator  F=(£p+£)  Itself  in  an  analogous  way  would  then  yield 
the  equations  for  separating  the  elastic  and  inelastic  effects  as  first 
order  differential  equations.  The  operator  F with  such  properties,  which  can 
be  called  an  adiabatic  momentum  operator,  can  be  formally  defined  by 


(3.11) 

^(P^)  = 0.  n=0, 1 

Applying  this  F twice  to  j,  yields  the  same  result  as  applying  the  F*  given 
by  (3.9).  From  (3.11)  one  finds 


y+f&+  + ifpa"  + p(y+g+  + y'a-)  = o (3.12) 

The  above  constraint  is  clearly  not  the  usual  variation  of  constant  con- 
straint *'“,s  (namely  fhat  U+  P ft4  + U pg,  equals  zero). 

Using  (3.1 2)  and  (3.2),  Eq.  (3.1)  can  be  reduced  to  a set  of  first 
order  differential  equations  for  the  a*; 
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(3.13) 


Since  the  coupled  equations  (3.13)  have  real  coefficients  a purely  real 
set  of  solutions  can  be  constructed.  Since  a solution  post-multiplied  by  a 
constant  matrix  is  also  a solution  of  (3.13)  the  two  real  matrix  solutions 
can  be  used  to  generate  a solution  with  any  desired  boundary  conditions, 
including  the  complex  ones.  Eq.  (3.13)  can  then  be  reduced  to  an  equation  of 
the  same  size  with  real  coefficients. 


- - 

a • • 0 

- - 

a 

BPA  - BPA  BPB  - BPB 

a 

d 

dR 

- i 

(3.14) 

b 

APA  - APA  APB  - APB 

b 

where 


U±  = (l/2)'/j(|±  IA)  a~  « (1/2)t/?  (b  + la) 


(3.15) 


and  the  matrix  P (purely  imaginary  for  a real  basis  set  X)  is  defined  in  (2.4). 
The  above  equations  are  convenient  for  numerical  work  since  their  Inte- 
gration requires  real  arithmetic  only. 
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IV.  A NUMERICAL  EXAMPLE  - TRANSITION  PROBABILITIES 
A&D  THE  SCATTERING  VAVEFUNCTION  FOR  AN  ATOM-ION  COLLISION 


In  this  section  a numerical  solution  of  the  Eq.  (3.14)  derived 
In  Sec.  II  is  presented.  The  transition  probabilities  and  the 
scattering  wavefunction  are  evaluated  for  an  atom-ion  col  I is  ion 
In  the  two-state  approximation. 

The  physical  system  chosen  is  the  excitation  collision 

He+  + Ne(2p6)  - He+  + Ne#(2p53s)  (4.1) 


at  collision  energy  of  70.9  e. v.  (2,60566  Hartrees).  The  above  inelastic 
process  has  been  studied  extensively  In  the  past  decade,  both  computat ional - 
lyiS'ts.tt'  and  exper  i mental  ly.  **“  31  Thema  in  features  of  the  experimental 
results  can  be  explained  by  a two-state  curve  crossing  model  in  which  the 

coupling  results  from  the  radial  part  of  the  Hamiltonian  only,  so  that  the 

ou 

angular  momentum  Is  preserved  throught  the  collision.  The  coupled  equations 
In  the  diabetic  representat ion  which  describe  mathematically  the  process 
with  the  above  approximations  can  be  written  as  19 


(<•  2) 


The  potential  matrix  elements  used  are  those  of  ref.  15. 
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Eq.  (4.2)  can  be  converted  Into  two  coupled  equations  in  the 
adiabatic  representation  by  solving  the  algebraic  problem  of  diagonal iza- 
tlon  of  the  matrix  ¥,  which  for  two  states  can  be  performed  analytically. 
The  transformation  matrix  is  given  by  (A. 7)  and  the  resulting  set  of 
adiabatic  equations  can  be  written  in  matrix  form  as 


where 

W.,2  - J (V.,  + V12)  + [(  V 1 ,-V  1 a ) */  4 + V 1 2 *]  ■/«  (4.4) 

and 

P„  = l(dt/dR)'/2/(l+t*>,  t = ( V22-V, ,)/ 2Vj2  (4.5) 


The  method  from  Sec.  2 can  now  be  applied  to  (4.3)  directly, 
corresponding  uncoupled  adiabatic  equations  are  given  by 


r 

0“ 

+ 

w 11 0 

0 1 1 

0 12 

- E 

0 1 1 

0 12 

L° 

p2_ 

.0  W2Z. 

_0  2 1 

02  2. 

_02  1 

02  2_ 

The 


(4.6) 


The  solutions  of  (4.6)  which  are  needed  for  solving  (4.3)  were  constructed 
numerically  using  Gordon's  method. *t52  Any  desired  solution  of  (4.6)  is 
then  written  as 
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• 

1 1 1 

A| 

O" 

*1 1 *12 

B| 

o' 

bM 

b.f 

tn 

♦ 22 

0 

a z 1 *22 

+ 

0 

82 

b2  1 

. b22 

where  Aj  and 

Bj  are  two 

linearly  independent  real 

solutions  of 

each 

(4.7  ) 


(4.6):  A.  is  the  solution  which  decays  exponentially  in  the  classically 
forbidden  region  and  is  given  by  A*f(  l/2kj)'/2  cos  (k.r  asympot  ical  1 y 

in  the  classically  allowed  region.  B.  is  the  other  linearly  independent 
solution  normalized  to  a unit  Wronsklan,  which  is  given  asympot ical ly  in 
the  classically  allowed  region  as  Ba/(  l/2k.)'^2  sin  (k{r  - 5^. 

The  initial  conditions  on  the  matrices  a and  b used  in  the  calculations  are 


a = 


ai 1 ai 2 

1 o' 

b|l  bl2 

o’ 

.a 

= 

b = 

= 

a 2 1 8 2_2 

0 !_ 

b2l  b2  2. 

9 oj 

(4.8) 


The  integration  of  the  equations  was  started  at  a value  of  the  radial  coor- 
dinate low  enough  to  obtain  converging  results  for  the  transition  probabilities 
to  three  significant  figures.  Starting  the  inelastic  integration  at  0.1  a.  u.  to 
the  left  of  the  outermost  turning  point  was  usually  sufficient.  The  integra- 
tion was  stopped  at  a value  of  R = Rf  where  the  matrices  a and  b became 
constant  to  the  same  accuracy.  (-*5  a.u.)  The  matrices  a+(Rf)  and  a”(R^)  were 
then  constructed  from  the  final  values  of  the  matrices  a and  b according 
to  (3. 15)  and  an  inelastic  scattering  matrix  S'  was  constructed 

a*) 

The  conventional  scattering  matrix  S can  be  obtained  in  terms  of  the  matrix 
S' by  comparing  their  definitions  (2.1j)  and  (4.9).  since  the  reference 
solutions  Ffj(2.Jl»)  and  U^-,  ( 3.5)  differ  by  the  elastic  phase  shift 
the  relation  between  S and  S'  Is  given  by 

(4.10) 


i 
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where  T|  Is  a diagonal  matrix  with  elements 

T){  f = exp  ( lTlj)  (4.11) 

- The  numerical  integration  was  performed  by  solving  (3.14  ) using  a 
standard  coupled  differential  equations  Integration  subroutine.32  The 
results  of  the  calculation  and  several  angular  momenta  for  total  energy 
of  the  collision  of  70.9  e.  v.  (2.60566  A.  U. ) are  given  in  the  second  column 
of  table  1.  The  third  column  numbers  are  from  Ref.  28  and  have  been  obtained 

by  solving  the  equivalent  set  of  equations  in  the  diabatic  representation  (4.2) 

The  elements  of  the  first  row  of  the  matrices  a and  b as  a function  of  R 

for  the  actual  imposed  initial  condition  in  the  numerical  Integration, 

(Eq.  (4.8))  are  presented  graphically  on  part  b of  Figure  1.  Those  are  the 
quantities  calculated  as  a function  of  R in  the  actual  numerical  integration  and  the 
virtual  absence  of  any  oscillations  or  instabilities  is  evident  from  the 
figure. 

The  first  row  of  the  wavefunctlon  on  the  l.h.s.  of  (4.6)  as  a function 
of  R in  the  vicinity  of  the  "crossing  point"  is  presented  graphically  on 
part  c of  the  Figure.  The  highly  oscillating  wave  functions  shown  here 
are  the  quantities  which  would  have  to  be  foPowed  numerically  by  the 
integration  routine  if  a solution  of  (4.3)  had  been  attempted  directly 
for  the  same  initial  condition. 


IB 


The  physically  meaningful  boundary  condition  on  the  matrices  g+  and 
g"  at  large  R Is  that  the  matrix  g"  ( R^)  be  a unit  matrix.  In  that  case, 

a 

the  first  and  second  rows  of  the  matrices  rj+  and  g"  for  every  R correspond 
to  the  solution  for  which  the  particle  is  incoming  in  the  lower  or  upper 
channel,  respectively.  The  matrices  g+  and  g"  corresponding  to  such 
boundary  conditions  which  will  be  denoted  with  c+ .and  c-can  be  obtained  by 
post-multiplying  the  matrices  g*  and  g"  for  every  R with  the  Inverse  of 
the  final  value  of  the  matrix  g“  namely 

(4. 12) 

The  elements  of  the  first  row  of  the  wavef unction  ji  corresponding  to  such 
boundary  condition  can  be  obtained  at  every  R by  pre-mul tipi ying  the  matrices 
c+  and  c“  with  the  reference  solution  matrices  U4  and  U~  as  defined  by  (3.5), 
respect  ivel y, and  are  shown  on  part  d of  the  figure. 


I 
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V.  DISCUSSION 

• ] 

The  method  for  treating  inelastic  collisions  presented  here  Is  based 

on  an  explicit  separation  of  the  elastic  and  inelastic  effects  of  the  col- 
lision. It  is  shown  that  the  oscillatory  structure  of  the  coupled  equations, 
j which  is  one  of  the  main  sources  of  dif f icul t ies • in  developing  a solution 
method,  is  almost  entirely  due  to  elastic  effects.  The  solutions  of  the 

elastic  part  of  the  problem,  which  can  be  obtained  relatively  easily,  can  j 

then  be  used  to  derive  a set  of  equations  for  the  purely  inelastic  effects. 

The  variables  in  the  latter  are  a measure  of  the  probability  of  finding  the  | 

particles  in  a given  internal  state  at  any  point  in  the  collision  and  the  j 

virtual  absence  of  oscillations  in  their  change  with  the  interparticle  dis- 
tance in  the  model  problem  treated  here  is  evident  from  the  numerical  cal- 
culations presented. 

The  present  formalism  can  be  shown  to  be  the  exact  analog  of  some 
semlclassical 3 * and  uniform  semicl ass  leal  27  treatments  in  the  literature. 

To  obtain  the  latter  one  only  needs  to  substitute  the  semiclassical  or 
uniform  semiclassical  approximation  to  the  elastic  wavef unct ion:  the  equa- 
tions for  the  Inelastic  effects  remain  the  same.  In  such  a way  the  present 
formalism  offers  a simple  way  of  deriving  the  above  approximations  and 
also  allows  one  to  express  the  final  equations  in  simple 
form. 

Introducing  some  additional  classical  limit  approximations,  as  shown 
In  ref.  3*,  reduces  the 'present  equations  to  the  classical  path  equations, 
which  forms  a convenient  framework  for  intuitive  interpretation  of  the  variables 
In  the  Inelastic  equations  In  terms  of  their  classical  limit  behavior.  A 
more  detailed  discussion  of, the  semiclassical  aspects  of  the  problem  Is 
given  elsewhere.  *•*  * 
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APPENDIX. 

I.  Transformations  between  the  equations  in  different  representations 
The  coupled  equations  (2.8,  2.9)*  rewritten  here 

(Ip  + P)2jj,  + Vi£  ■ Ejj>  (A.  I) 

can  be  converted  into  the  adiabatic  representation  using  the  unitary 
transformation  M which  diagonalizes  the  potential  energy  matrix  V 

MVM  - Va  (A. 2) 

to  transform  the  solution  matrix  <£ 

f (A. 3) 

The  tilde  symbol  ' above  denotes  a hermition  conjugate  matrix. 

The  transformed  equations  in  the  adiabatic  representation  are  obtained  by 
substituting  (A. 2)  and  (A. 3)  in  (A.l)  and  premultiplying  the  resulting 
equation  by  M,  yielding 

(Ip  + Pa)j&a  + VY  - E^a  (A.4) 

where 

P*  - fiCjp  + £)H  (A. 5) 

Similarly,  the  coupled  equations  (A.l)  can  be  converted  into  the  diabatic 
representation  by  means  of  the  unitary  transformation  which  diagonalizes  the 
kinetic  energy  part  of  the  Hamiltonian 

N(lp  + P)2N  - Ip*  (A. 6) 

A solution  of  the  second  order  differential  equation  (A. 6)  is  given  by  the 


solution  of  the  first  order  differential  equation 
PN  - £N 


(A.  7) 
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Eq.  (A. 7)  gives  a rigorous  definition  of  the  diabetic  representation 
for  the  case  of  atomic  collision  problems  where  It  cannot  be  easily 
defined* In  terms  of  the  asymptotic  basis  sets.5  However,  since  the  trans- 
formation or  a set  of  equations  In  the  adiabatic  representation  Into  the 
diabatlc  representation  requires  solving  a coupled  set  of  differential 
equations  (A.7),  It  Is  practically  never  used. 

The  definition  of  the  adiabatic  and  the  diabatlc  representation  given 
in  Sec.  2 is  based  on  the  choice  of  the  internal-orientational  basis  set 
that  leads  to  a particular  representation.  Alternatively,  the  adiabatic 
representation  has  been  defined3  by  the  vanishing  of  the  off-dlagonal  elements 
of  the  matrix  V,  and  the  diabatlc  one  by  the  vanishing  of  the  diagonal 
elements  of  P in  (2.  9).  when  the  basis  set  X is  complete  the  two  alternative 
definitions  are  ident ical.  5* 35  When  the  basis  set  used  to  derive  the  equa- 
tions Is  not  complete  they  are  however  different.  The  "adiabatic"  set  of 
equations  obtained  by  diagonalizing  the  potential  matrix  in  a finite  diabatic 
set  of  equations  corresponds  to  using  an  approximate  expansion  of  the  eigen- 
functions of  the  internal  part  of  the  hamiltonian  at  every  R In  terms  of 
the  asymptotic  internal  states.  The  "diabatic"  set  of  equations  obtained 
by  diagonalizing  the  kinetic  energy  matrix  operator  in  a finite  adiabatic 
set  of  equations  similarly  corresponds  to  equations  obtained  by  using  an 
approximation  to  the  asymptotic  internal  wavefunct ions  in  terms  of  the  local 
eigenfunctions  of  the  internal  hamiltonian  at  every  R,  and  it  is  weakly 
dependent  on  R. 

2.  Symmetry  Properties  of  the  Coupled  Equations 

Eq.  (2.  B ) Is  a self-adjoint  second  order  matrix  differential  equa- 
tion In  the  sense  that  the  hermitian  conjugate  of  the  solution  Is  an  Inte- 
grating factor  for  the  equation  and  can  be  used  to  reduce  its  order  by  one.  **• 


22 


Premultiplying  (2.11)  by  the  hermitian  conjugate  of  the  solution,  taking  the 
hermitian  conjugate  of  the  matrices  on  both  sides,  end  subtracting  one  gets 

“ (p*|)|  + 2rjppJ  + (p|)PJ  +1(pP)|)  = 0 (a. 8) 

where  the  explicit  form  of  T (2.9)  has  been  used. 

Rearranging,  (A. 8)  can  be  written  as 

P[fpJ  - (p|)|  + 2(|p£)]  = PF('|,|)  = 0 (A-9) 

or  integrating  both  sides,  after  substituting  p - -i(d/de) 

1 P | - (pjpi  + 2(|>J)  = F(]f,|)  = const  (A.  10) 

Hence,  the  hermitian  conjugate  of  the  solution  of  (2.11)  Is  an  integrating 
factor  for  the  same  equation  and  the  equation  is  self-adjoint.  Equation 
(A.  10)  can  be  viewed  as  the  analog  of  the  flux  conservation  relation  for 
the  Schrodinger  equation.  It  can  be  obtained  from  the  flux  conservation 
relation  for  the  Schrodinger  equation  simply  by  substituting  the  basis 
set  expansion  for  the  wavef unction  and  integrating  over  the  r,Q  dependence.  5 
The  value  of  the  constant  matrix  F can  be  evaluated  at  R=0  where  the  wave- 
function  matrix  $ vanishes  and  it  is  zero. 

The  self-adjoint  property  of  the  coupled  equations  or  their  "Hermicity" 
has  been  a point  of  some  confusion  in  the  literature.  It  was  argued  (c. f. 
ref.  3,  sec.  9)  that  the  matrix  B (2.7)  is  not  hermitian  and  therefore 
the  coupled  equations  are  not  Hermitian.  While  the  matrix  B does  indeed 
represent  a non  self-adjoint  term  In  (2.  9)  the  matrix  operator  Pp  Is  also 
not  self  adjoint.  The  easiest  way  to  recognize  this  point  Is  to  notice  that  Pp 
Is  a product  of  two  none ommut ing  self  adjoint  matrix  operators  so  that  the 
whole  Hamiltonian  matrix  operator  is  self-adjoint,  as  one  would  expect. 
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Table  I: 

The  Transition  Probabilities  |S|2|*  for 

He  + Ne  _ He  + Ne#at  collision 

• 

energy  of  70.9  ev  and  different  angular 

momenta 

Angular  Momentum  |Siz|2(this  work) 

|S|||*(ref.  2 B) 

200 

0.215 

0.215 

300 

0. 146 

0. 146 

310 

0.202 

0.202 

320 

0.0941 

0.  094 1 

330 

0.0226 

0.0226 

Legend  to  Figure  I 


The  adiabatic  potential  energy  matrix  elements  W end  the  coupling  matrix 


element  P(2  ere  shown  in  (a)  as  a function  of  the  radiat  coordinate  R. 

The  W | | , WI2  and  E are  in  Hartrees,  while  P||/10  is  in  bohr-1.  The  first  row 
^elements  of  the  matrices  a,  solid  line,  and  b,  dashed  line,  obtained 
by  solving  numerically  (3.14)  with  the  initial  conditions  (4.6)  are  shown  in 

(b) .  The  wavefunction  (4.5)  for  the  same  initial  conditions  is  shown  in 

(c) .  The  actual  wavefunction,  corresponding  to  no  incoming  wave  in  the 
upper  channel  is  shown  in  (d),  where  the  solid  line  describes  the  real 
part  and  the  dotted  line  the  imaginary  part  of  the  wavefunction.  In  each 
case  (b)  to  (d),  the  upper  curve  or  curves  refer  to  channel  2 (the  more 
highly  internally  excited  channel),  and  the  lower  curve  or  curves  refer 
to  channel  1.  The  structure  of  the  adiabatic  potential  matrix  elements 
around  the  avoided  crossing  is  shown  in  the  insert  on  part  a. 


i 


'f  M - : #y»S  - 


472 :GAK : 716 : tam 
78u472-608 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  GEN 


• No. 

Copies 

y 

Office  of  Naval  Research 

600  North  Quincy  Street 

Arlington,  Virginia  22217 

Attn:  Code  472  2 

j 

ONR  Branch  Office 
536  S.  Clark  Street 
Chicago,  Illinois  60605 
; Attn:  Dr.  George  Sandoz  1 

ONR  Branch  Office 
715  Broadway 

New  York,  New  York  10003 

Attn:  Scientific  Dept.  1 

ONR  Branch  Office 
1030  East  Green  Street 
„ Pasadena,  California  91106 

Attn:  Dr.  R.  J.  Marcus  1 

r* 

ONR  Area  Office 

One  Hallidie  Plaza,  Suite  601 

San  Francisco,  California  94102 

Attn:  Dr.  P.  A.  Miller  1 

» 

ONR  Branch  Office 

Building  114,  Section  D 

666  Summer  Street 

Boston,  Massachusetts  02210 

Attn:  Dr.  L.  H.  Peebles  1 

Director,  Naval  Research  Laboratory 

Washington,  D.C.  20390 

Attn:  Code  6100  1 

The  Assistant  Secretary 
of  the  Navy  (R,E&S) 

Department  of  the  Navy 
Room  4E736,  Pentagon 

Washington,  D.C.  20350  1 

Commander,  Naval  Air  Systems  Coomand 

Department  of  the  Navy 

Washington,  D.C.  20360 

Attn:  Code  3 10C  (H.  Rosenwasser)  1 


No. 

Copies 

Defense  Documentation  Center 

Building  5,  Cameron  Station 

Alexandria,  Virginia  22314  12 

U.S.  Army  Research  Office 
P.0.  Box  1211 

Research  Triangle  Park,  N.C.  27709 
Attn:  CRD-AA-IP  1 

Naval  Ocean  Systems  Center 

San  Diego,  California  92152 

Attn:  Mr.  Joe  McCartney  1 

Naval  Weapons  Center 

China  Lake,  California  93555 

Attn:  Dr.  A.  B.  Amster 

Chemistry  Division  1 

Naval  Civil  Engineering  Laboratory 

Port  Hueneme,  California  93401 

Attn:  Dr.  R.  W.  Drisko  1 

Professor  K.  E.  Woehler 
Department  of  Physic**  & Chemistry 
Naval  Postgraduate  School 
Monterey,  California  93940  1 

Dr.  A.  L.  Slafkosky 
Scientific  Advisor 
Connandant  of  the  Marine  Corps 
(Code  RD-1) 

Washington,  D.C.  20380  1 

Office  of  Naval  Research  * 

800  N.  Quincy  Street 

Arlington,  Virginia  22217 

Attn:  Dr.  Richard  S.  Miller  1 

Naval  Ship  Research  and  Development 
Center 

Annapolis,  Maryland  21401 
Attn:  Dr.  C.  Bosmajian 

Applied  Chemistry  Division  1 

Naval  Ocean  Systems  Center 
San  Diago,  California  91232 
Attn:  Dr.  S.  Yamamoto,  Marine 

Sciences  Division  1 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  359 


No. 

Copies 


472  :GAN : 716 : tan. 
78u472-608 


No. 

Copies 


Dr.  Paul  Delahay 
New  York  University 
Department  of  Chemistry 
New  York,  New  York  10003 

Dr.  R.  A.  Ostcryoung 
Colorado  State  Univeraity 
Department  of  Chemistry 
Fort  Collins,  Colorado  80521 

Dr.  E.  Yeager 

Case  Western  Reserve  University 
Department  of  Chemistry 
Cleveland,  Ohio  41106 

Dr.  D.  N.  Bennion 
University  of  California 
Chemical  Engineering  Department 
Los  Angeles,  California  90024 

Dr.  R.  Av  Mar/41  a 

CaliforniaVPhstitute  of  Technology 
Dcpartmeny&f  Chemistry 
Pasadena,'  California  91125 

Dr.  J.  J.  Auborn 

Bell  Laboratories 

Murray  Hill,  New  Jersey  07974 

Dr.  Adam  Heller 

Bell  Telephone  Laboratoriea 

Murrary  Hill,  New  Jersey  07974 

Dr.  T.  Katan 

Lockheed  Missiles  6 Space 
Co , Inc . 

P.0.  Box  504 

Sunnyvale,  California  94088 

Dr.  Joseph  Singer,  Code  302*1 
NASA- Lewis 

21000  Brookpark  Road 
Cleveland,  Ohio  44135 

Dr.  B.  Brumaer 
BIC  Incorporated 
Five  Lee  Street 

Cambridge,  Massachusetts  02139 


Library 

P.  R.  Mallory  and  Company,  Inc. 
Northwest  Industrial  Park 
Burlington,  Massachusetts  01803 

Dr.  P.  J.  Hendra 
University  of  Southhampton 
Department  of  Chemistry 
Southhampton  S09  5NH 
United  Kingdom 

Dr.  Sam  Perone 
Purdue  University 
Department  of  Chemistry 
West  Lafayette,  Indiana  47907 

Dr.  Royce  W.  Murray 
University  of  North  Carolina 
Department  of  Chemistry 
Chapel  Hill,  North  Carolina  27514 

Naval  Ocean  Systems  Center 
San  Diego,  California  92152 
Attn:  Technical  Library 

Dr.  J.  H.  Ambrua 
The  Electrochemistry  Branch 
Materials  Division,  Research 
4 Technology  Department 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  Maryland  20910 

Dr.  6.  Goodman  • 

Globe-Union  Incorporated 
5757  North  Creen  Bay  Avenue 
Milwaukee,  Wisconsin  53201 

Dr.  J.  Boachlar 
Electrochimica  Corporation 
Attention:  Technical  Library 
2485  Charleston  Road 
Mountain  View,  California  94040 

Dr.  P.  P.  Schmidt 
Oakland  University 
Department  of  Chemistry 
Rochester,  Michigan  48063 


1 


1 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  359 
No. 

Copies 

Dr.  H.  Richtol 

Chemistry  Department 

Rensselaer  Polytechnic  Institute 

Troy.  New  York  12181  1 

Dr.  A.  B.  Ellis 

Chemistry  Department 

University  of  Wisconsin 

Madison,  Wisconsin  53706  1 

Dr.  M.  Wrighton 
Chemistry  Department 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139  1 

Larry  E.  Plew 

Naval  Weapons  Support  Center 

Code  3073,  Building  2906 

Crane,  Indiana  47522  1 

S.  Ruby 
DOE  (STOR) 

600  E Street 

Washington,  D.C.  20545  1 

Dr.  Aaron  Wold 

Brown  University 

Department  of  Chemistry 

Providence,  Rhode  Island  02192  1 

Dr.  R.  C.  Chudacek 
McCraw-Edison  Company 
Edison  Battery  Division 
Post  Office  Box  28 

Bloomfield,  New  Jersey  07003  1 

Dr.  A.  J.  Bard 

University  of  Texas 

Department  of  Chemistry 

Austin,  Texas  76712  ^ 

Dr.  M.  M.  Nicholson 
Electronics  Research  Center 
Rockwell  International 
3370  Mlraloma  Avenue 

Anaheim,  California  92803  1 

Dr.  M.  0.  Sceats 
University  of  Rochester  . 

Department  of  Chemistry 

Rochester,  New  York  14627  1 


472 :GAN: 716: tar 
78u472-608 


